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Abstract. We introduce a special class of random matrices (DUE) whose 
spectral statistics corresponds to statistics of microscopical quantities detected 
in vehicular flows. Comparing the level spacing distribution (for ordered 
eigenvalues in unfolded spectra of DUE matrices) with the time-clearance 
distribution extracted from various areas of the flux-density diagram (evaluated 
from original traffic data measured on Czech expressways with high occupancies) 
we demonstrate that the set of classical systems showing an universality associated 
with Random Matrix Ensembles can be extended by traffic systems. 
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Connection between Random Matrix Ensembles and certain transport systems is not new (see e.g.[l] for 
a general overview). Indeed, in Ref.[2, 3] authors drawn the attention to the bus transportation system 
in Cuernavaca, Mexico, where a peculiar transfer of information inside the system led, surprisingly, to 
the universal configuration of gaps among buses. The rigorous study [4] confirmed that the detected 
link between Cuernavaca buses and Gaussian unitary ensemble (CUE) is not accidental. Authors of the 
research formulated a relevant microscopic model-scheme and analytically proved that the time headway 
distribution of buses conforms to the level spacing distribution for CUE. This surprising knowledge 
attracted an attention of scientists (e.g. [1]) since systems connected to Random Matrix Theory are 
universal (in a certain sense). 

Until now, natural endeavor to interconnect Random Matrix Theory (RMT) with vehicular systems 
has not led to a success. However, a partial progress has been achieved in [5, 6, 7, 8, 9] where it is 
proved that microscopical arrangement of vehicles can be predicted by means of a certain one-dimensional 
gas inspired by the Dyson’s gases that are well-known in RMT. Moreover, another attempts to describe 
(analytically) a microstructure of vehicular ensembles with help of statistical instruments (e.g. [10, 11, 12]) 
led us to believe that our effort may be successful. Therefore, the main goal of this work is to find (and 
analyze) a new class of random matrices whose spectral properties correspond to a micro-structure of 
real-road traffic samples. 

This paper is organized as follows. In the first section we introduce and analyze the class of the 
so-called damped matrices (DUE). We are focused on the level density, the procedure of the unfolding, 
and the level spacing distribution. Connection to the theory of one-dimensional traffic gases is discussed 
in the second section. Section Three brings in a comparison between the level spacing distribution of DUE 
matrices and the time-clearance distributions extracted from expressway data samples. This is followed 
by corresponding discussion and conclusions in the final section. Subsidiary derivations are included in 
the mathematical Appendix. 

1. DUE Damped unitary ensembles 

Being inspired by the work [13] (studying Calogero-Moser models with various potentials) we introduce 
the (/—parameterized class DUEg(A^) whose N x N matrices Hg = {hkj{g ))^fulfil the following axioms: 

(i) the parameter g ^ 0 is fixed; 

(ii) elements hkk{g) are chosen as independent Gaussian variables with zero mean 
and unit variance, i.e. hii{g), h 22 {g), • ■ •, hNN{g) ^ N(0,1) are i.i.d.; 

(iii) if fc ^ j then 

. / ^ __ 27rig _ 

^ A^sinh[(27r(fc — j)/iV] ’ 

i.e. the off-diagonal elements are deterministic. 

Owing to the definition the off-diagonal elements are purely imaginary complex numbers whose 

absolute value become smaller with distance from the diagonal. Therefore we refer such a class to as the 
Damped Unitary Ensembles. In the context of the article [13] the matrices in DUEg(iV) represent simplified 
versions of the N x N Lax matrices derived for hyperbolic Calogero-Moser model (see [13] for details). 

1.1. Level density for DUE matrices 

Let <j{Hg) be the spectrum of Hg, i.e. cr{Hg) = S C : /c = 1, 2,..., N} is a set of ordered eigenvalues 
X(i) ^ X( 2 ) ^ ... ^ X( 7 v). Since Hg = H^, where H^ is Hermitian conjugated matrix, one finds that 
(j{Hg) C R and therefore the above-referred ordering is correctly defined. Now we can introduce the 
probability density Pkix(k)) for fcth eigenvalue x^k) and the level density 

1 ^ 

W7v(x) = 


( 1 ) 
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Figure 1. Graphs of the level density Histograms have been determined 

for five thousand realizations of DUE(256) matrices with parameters g = 0.2 (a 
northwest sub-figure), g = 0.5 (a northeast sub-figure), g = 1.5 (a southwest sub¬ 
figure), and g = 2.5 (a southeast sub-figure). Solid curve (blue)/dashed curve 
(magenta) represent the approximation (2) for values e, 'd estimated with help of 
CSC/KE methods, respectively. 


Unfortunately, the famous Semi-Circle Law (derived in [14] for the level density of classical random matrix 
ensembles) is not applicable for DUE matrices. Instead of a circular shape the level density of damped 
matrices conforms to a probability density taken from two-parametric family 




m [ exp[-j|^] ... |a:|<£, 

e \ 0 ... |x| > e, 


( 2 ) 


where C, exp[—dx ensures the proper normalization. Indeed, numerical tests show (as is 
illustrated in figure 1) that the level density for large N is very accurately approximated by the function 
qs,e{x) for estimated values e = e{g, N),!} = ^^{g, N) of the so-called critical eigenvalue e and curvature i}. 
Both the presented estimation-procedures are based on principles of MDE (minimum distance estimation). 
The first of them (CSE - chi-square estimator) minimizes the statistical distance 


xA,£) = ^J IHA - (3) 

where h{x) is an empirical histogram. For potential reconstructive purposes we are specifying that the 
set S consisting of 1 280 000 eigenvalues (taken from 5 000 realizations of 256 x 256 matrices) generates 
the domain Dom(/i) = [min(S'), max(5')] of the empirical histogram. In our numerical tests, this domain 
Dom(/i) has been divided into 200 equidistant sub-intervals. 

The second statistical test is based on the standard Kolmogorov estimator (KE) minimizing the 
supremum of the absolute difference between the estimated distribution function and the empirical 
distribution function 

1 ^ 

Hn{x) =^x). (4) 

^ k=l 


Here, I{x) stands for the indicator function. The optimal values of estimated parameters are compared 
graphically in figure 2. Marked discrepancies near the origin can be attributed to the following facts. If 

g = 0 then one can trivially express the level density as 2 /v27r and our estimation procedure tries, in 
fact, to approximate the normal distribution by a function q^^^(x). Moreover, 


Vd > 0 : lim ^(x) 

6 —^0-j- 


lim _ e 27^ 

^-^ 0 + \r2/K>t 


6{x), 
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where S(x) € is the Dirac function (see Appendix 5.1). Since both these distributions tend (for small 
values of parameters) to the same generalized limit it is difficult to distinguish them numerically. It means 
that estimations used (applied for small values of g') are highly sensitive on a type of estimators, quality of 
analyzed data and so on. On the other hand, nuances in parameters estimated by different methods cause 
only imperceptible deviations in curves of the level density (see solid and dashed lines in the northwest 
sub-figure in figure 1)). 



Figure 2. Estimates of the critical eigenvalue e and curvature '& in (2). A curve 
and dashed line correspond to i. and ^ obtained by CSE method. Squares and 
circles display the estimates e and ^ calculated by means of KE. 


1.2. Procedure of the unfolding 


Before statistical analysis of random-matrix spectra the standard methodology of Random Matrix Theory 
requires a transition to the so-called unfolded spectrum. The main benefit of that transition is an unified 
view into a structure of heterogeneous systems. To be specific, the procedure of the unfolding converts 
eigenvalue-spectrum cf{Hg) with an arbitrary level density ujn{x) (depending on specific properties of a 
system) to the unfolded spectrum a{Hg) = {i(fe) £ C : k = 1,2,..., N} in which unfolded eigenvalues are 
uniformly distributed in (0,iV). Thus, the level density of the unfolded spectrum is associated with the 
uniform distribution U{0,N). Mathematically, the unfolding is a procedure consisting of two operations: 
an individual mapping i—>■ defined by 

r^(k) 

X{k) := / WAr(?/)dy, (5) 


and a global scaling 


^{k) • — 


X(k){N - 1 ) 

Efc=i (^(fc+i) “ ^(fc)) 


( 6 ) 


According to the theorem on probability integral transformation (see [15]) it holds x ^ U(0,1), which 
results (for large N) in the fact x ^ \J{0, N). Moreover, the average gap between two succeeding unfolded 
eigenvalues is equal to one. 

Specifications: All the previous considerations require a large value of N. Thus, in our numerical 
realizations we consistently use N = 256. 


1.3. Level spacing distribution 

In this subsection we aim to deduce an analytical formula predicting the probability density for gaps 
Sk ■= i(fe+i) — i(/c) between the ordered eigenvalues ^ X( 2 ) ^ ... ^ X( 7 v) in the unfolded spectrum 
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a{Hg). Here, such a function is referred as the level spacing distribution (LS-distribution) and denoted by 

p(s)- 

The most trivial variant of ensembles investigated is the class DUEo(-/V) of totally damped matrices. 
The matrices H G DUEo(iV) are diagonal and therefore the corresponding level density reads ujn{x) = 
_^Q-x /2 jy Thus, the first stage of the unfolding-procedure 

^ dy = - (^1 + Erf (7) 

produces i.i.d.-numbers from the uniform distribution 17(0,1). Hence, one can easily verify that level 
spacing distribution in DUEo(iV) can be expressed as 


p(s) = 0(s)e (8) 

In the article [13] authors derived that the joint probability density for eigenvalues in the hyperbolic 
Calogero-Moser model corresponds to the law 

N N 

P{xi,X2 ,..., xn) oc exp[—H /Co 

k=l k=l 

where 0,'!' are constants and ICm{x) stands for the Macdonald’s function of the order m - solution of 
the modified Bessel’s differential equation of the second kind x^y" + xy' — -I- m?)y = 0. Having regard 

to the fact that DUEg(iV) matrices can be understood as numerical simplifications of Lax matrices in the 
hyperbolic Calogero-Moser model (as is discussed in [13]), the function (9) can be considered as a suitable 
approximation for the joint probability density in the ensemble DUEg(iV). Although the expression (9) is 
difficult to handle, it allows to obtain a Wigner-type surmise for the level spacing distribution in DUEg(A^). 
Considering N = 2 (similarly to the original Wigner idea) we have 



1-f 


Airlg 


N{xj - Xk) 


,( 9 ) 


P{xi,X 2 ) OC e 


— 1^2 
^0 




(x2 - a:i)2y ■ 


( 10 ) 


From the general relation £{r) = <d{r) (^J^P{r + x,x) dx + J^P{x,r + x) dx) we can obtain the 
probability density 


e{r) oc 0(r)/Co 1 -|- ^ e (11) 

for a distance between non-unfolded eigenvalues in DUEg(2). Application of the approximation K.q{z) « 
+ 0{1/z) is resulting in the formula 

47 r^g 

£{r) oc re ~, (12) 

describing a behavior of the distribution studied at small distances r. An immediate consequence of this 
result is evidence of repulsions between eigenvalues that are rapidly stronger than those detected in the 
standard matrix ensembles (i.e. GDE, GUE). 

Combining the two asymptotic predictions (8) and (12) we suggest (analogously to [13]) the level 
spacing distribution in the form 


p(s) = 

(13) 

where (according to Appendix 5.2) 


3 - 


D = ^ + a+ -^-, 

(14) 

a + 1 


/ n \ — 2 — 

A-^=2U^\ /Co+i(2v^). 

(15) 


On contrary to the pioneering work [13] our paper intends to investigate the level spacing distribution 
of DUE-matrices in more details. As the curse of the level density near the values x = ±£ is extremely 
steep (see, for example, the southeast sub-figure 1), the numerical representation of the unfolding (using 
an approximate formula (2) for the level density) is unstable in those regions. For this reason we analyze a 
truncated spectrum dm{Hg) = {x[k) '■ k = m + l,2,...,N — m} whose elements are re-scaled so that the 
mean spacing is equal to one. Then the set {sm+i, Sm+ 2 , ■ ■ ■, SAr-m-i} of empirical spacings is subjected 
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Level spacing s Level spacing s 


Figure 3. Level spacing distribution. The gray histograms correspond 
to empirical LS-distributions of DUEg(256) obtained for one thousand matrix- 
realizations and for g = 0.2 (a northwest sub-figure), g = 0.5 (a northeast sub¬ 
figure), g = 1.0 (a southwest sub-figure), and g = 2.0 (a southeast sub-figure). 
The continuous/dashed curves display the prognosis (13), in which the parameters 
have been estimated by means of MDE/MLE estimators, respectively. For all 
variants the value m = 50 has been used. 


to standard procedures estimating optimal values a, /3 in the functional premise (13). Here we present two 
standard estimators: MLE - maximum likelihood estimation and MDE - minimum distance estimation. 
For these purpose we introduce the log-likelihood function 

N—m — l A^—771—1 . A^ —771 —1 

C{a,l5) = {N-2m-l)\n{A)+a ^ ln(sfc)-/3 ^ - - D ^ Sfc(16) 

k—m+l k—m+1 fc—m+1 


and the weighted statistical distance 


x(a,/3) := d[h,p] 



\h{s) - p{s; {a, I3))\ ds, 


(17) 


TTS^ 

where h{s) is an empirical histogram for spacings. We remark that the weight factor wls) := 
has been chosen to eliminate the influence of long spacings and increase the influence of spacings that are 
close to the mean value. In addition, argmaxs^o^(s) = 1 and (i[0,1] = 1, which can be understood (if d 
is now understood as a general functional metric) as a benchmark. As graphically illustrated in figure 3 
LS-distributions of DUE—matrices is in full compliance with the analytical prediction (13). The optimal 
values of estimated parameters, i.e. (d,/3)MLE = argmax£(a,/3), (d,/3)MDE = argminy(a,/I) are then 
compared graphically in figure 4. 


2. Local thermodynamical traffic model 


As one of the main goals of this article is to find a connection between the above-referred matrix ensembles 
and freeway traffic samples, it is beneficial to show that the functional premise (13) can be understood 
as a steady-state distribution of particle headways in the thermal-like socio-physical traffic model with 
short-ranged repulsions between particles (originally introduced in [16]). Indeed, if considering circular 
ensemble of particles whose mutual distances (between subsequent particles) are denoted by si, S 2 , ■ • ■, 
and if understanding such an ensemble as thermodynamical (with a thermal parameter ^ 0) we are 
able to derive (following the approaches presented in [16]) that the headway distribution in the ensemble 
with a repulsive potential energy 


Uf,{si,S2,...,Sn) = E ( ~ -Mln(Sfe)) 

fc=l '' 


(m^^O) 


(18) 
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Figure 4. Estimated values of parameters in the functional premise (13). We 
visualize the parameters a, (3 estimated by the maximum likelihood estimation 
(bullets) and minimum distance estimation (crosses), respectively. Values of 
the statistical distance x{a.^f3) between the probability density (13) specified for 
estimated parameters and empirical histograms are plotted in the third sub-figure. 
Values of the log-likelihood corresponding to the probability density (13) 

specified for estimated parameters are plotted in bottom sub-figure. 


takes the form p(s) = As^'^e~~e~^^. Such a probability density fully corresponds to the generalized 
inverse Gaussian distribution distribution (13) and moreover, it fulfills all the substantial criteria for 
acceptability of analytical clearance distributions presented in the articles [8, 9]. Indeed, for all values 
of /i ^ 0 and all ^ > 0 the function p(s) generates the origin plateau, which means that it obeys the 
condition q G R”*" lims_,.o^ s“‘^p(s) = 0. Furthermore, the function p(s) is belonging to the class of 
balanced distributions (see [9]), which (together with the mathematical criteria (Tl—T5) discussed in 
[9]) completes all theoretical requirements for vehicular headway distributions. 

Since the presented thermodynamical traffic model represents a mathematically well-defined (and 
analytically soluble) alternative to the more realistic traffic model (see the discussion on the pages 6,7 in 
[8]) we have now a theoretical permission to use the above-referred statistics for a mathematical description 
of headway distributions in realistic vehicular streams. 
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3. Headway statistics in various regions of flux-density map 


In this section we try to answer the fundamental question of this work: ”Is there any link between fluid 
vehicular locations and eigenvalues in random matrices (similar to the link between Mexico buses and 
CUE—matrices presented in [2, 3, 4, 1])?” 

3.1. On the origin of traffic data used 

Vehicle-by-vehicle data analyzed in the paper have been recorded at the Expressway R1 (also called the 
Prague Ring) in Prague, the Czech Republic. For our measurements (conducted in cooperation with The 
Road and Motorway Directorate of the Czech Republic) we have chosen expressway-segments where traffic 
is strongly saturated and where therefore one can regularly detect synchronized flows. Having regard to 
intentions of this article we eliminate all free-flow states (where car-car interactions are negligible and 
therefore traffic behaves almost like poissonian system). For distinguishing congested-traffic states from 
free-flow states (see [18, 19]) we have applied the clustering algorithm K-means (see [17]). Here the chosen 
alternative of the method is based on minimization with respect to velocity of vehicles. 

3.2. The traffic alternative for the spectral unfolding 

Data record (adapted into a form suitable for further processing) has been divided into sub-samples 
Tj = {tju lj 2 ) ■ ■ ■, tjm} containing m consecutive time clearances (netto time intervals between succeeding 
cars passing a traffic detector located at the chosen lane), where m is the sampling size. For each sub¬ 
sample there were calculated the local flux Jj and local density Qj. The method used for those calculations 
is identical to the standard method described mathematically in [8] (page 4). Again, we emphasize that all 
sub-samples belong to a regime of congested streams when average speed of cars is dramatically reduced 
compared to the maximum allowable speed. 

In order to avoid an undesirable mixing between different traffic constellations we apply (similarly 
to the procedure presented for random matrices in the subsection 1.2) the so-called traffic unfolding. 
Here, such a procedure is composed of two components: a re-scaling and segmentation by a flux- 
density window. To be specific, for a certain window-size (Aj, Ag) we introduce a flux-density window 
^(^') >7) := [e, e + Ag] X [ J, J -I- Aj] that represents a small rectangular sub-region inside the fiux-density 
map. Then the re-scaling is understood as a transformation of the set Tj into the set {Tji,Tj 2 , •. • 
where Tjk = Furthermore, the segmentation 

7(£', J) = {j ■ Jj € [J, J -I- Aj] A Qj G [e, e -I- Ag]} (19) 

is understood as a procedure selecting sub-samples associated with a chosen flux-density window W{g, J). 
It means that the final sample of unfolded clearances looks like 

T{q,J) = {Tjk : j & I{g,J) A fc = 1,2,... ,m}. (20) 

For completeness, we add that in this study there is consistently considered m = 50, Aj = 400veh/h, 
and Ag = 5 veh/km. The total number of analyzed clearances is circa 45 millions (7 record-files gauged on 
7 different detectors between 2013-09-16 and 2014-05-29). Among 45 millions of measured clearances 14 
millions have been classified (with help of the K-means clustering algorithm) as congested (or transitional). 

3.3. Changes in headway statistics induced by varying fluxes and densities 

In this subsection we aim to analyze a clearance distribution extracted from unfolded data-sets (20) 
associated to a fixed flux-density window W{g,J). Withal we try to test a hypothesis that clearance 
distributions detected in regions of congested traffic belong to the two-parametric family of distributions 

p(r) = (21) 

If proving that hypothesis we find a direct link between inner structures of vehicular ensembles and spectra 
of damped random matrices. 

The testing mechanism is as follows. Varying a flux-density window we estimate a optimal parameters 
a, /3 ^ 0 for which the statistical distance 

/ TTT 2 o 

x{a,^) := d[p,q] = J J —|p(r; (a,/?)) - ( 7 (t) | dr (22) 
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Figure 5. Time clearance distributions in vehicular flows (analyzed for fast lanes 
of expressways). Histograms stand for empirical distributions analyzed inside the 
following flux-density windows: a northwest sub-figure: PF(48, 2300), northeast 
sub-figure: IT(54,1900), southwest sub-figure: IT(58,1300), and southeast sub¬ 
figure: IT(62,1300). The continuous/dashed curves visualize functions (21) 
specified for parameters estimated by MDE/MLE, respectively. The crosses 
correspond to a nearest distribution of eigenvalue-gaps, i.e. to a level spacing 
distribution whose statistical distance to traffic histogram is minimum. 
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Figure 6. Color representation of values /3 estimated by minimum distance 
estimation (23). We plot values /3mde in various flux-density windows W{(>,J). 
Squares symbolize locations of window-centroids (p + Agl2,J + Aj/2). The 
rectangle at the bottom right schematically visualizes the window-size used. Gray 
dots correspond to locations of individual data sub-samples before a 

separation procedure. Unstained areas indicate either a lack of data or belonging 
to a free-traffic regime. The plotted values are valid for a fast lane. 
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between a function (21) and clearance histogram g(r) is minimum. Practically, after eliminating outliers 
Tjfc > 6 and after a repeated re-scaling (so that the average clearance is equal to one again) we solve a 
minimization problem 


(d,/3)MDE = argmin„,/3 


re 


" |p(t; (a,/3)) - q(r)| dr, 


(23) 


where the histogram qlr) is enumerated for the data set T{q,J) = {fi : i = 1,2,..., M}, for which 
{T{g,J)) = 1. Analogously, the maximum likelihood estimations were carried out. For these aims we 
have defined an estimator 


By a maximization 


£(«,/?) 


M M 

M ln(A) + a ^ In(fi) - ^ 

i=l i=l 


1 

fi 


{a, /3)mle = 


argmaxa./3/l(a, /3) 


we obtain alternative values of estimated parameters. 


M 




(24) 


(25) 


log(a 


MDE' 



Traffic density p 


Figure 7. Color representation of values a estimated by minimum distance 
estimation (23). We plot values log(d]viDE) in various flux-density windows 
W{q^ J). Squares symbolize locations of window-centroids (^ + ^qI‘2^ J -h Aj/2). 
The plotted values are valid for a fast lane. 


In several figures we now depict selected results of our investigations. Method MDE provides 
estimations visualized in figures 6 and 7. Here one can recognize a behavior so typical for agent systems 
with intensive interactions among agents. Indeed, with increasing density a level of vehicular coordinations 
is growing up. This effect is clearly visible in figure 6 where estimated parameters /3 are continuously 
changing from negligible values (for regimes with small densities) to extremely high values detected for 
a strongly saturated traffic. Such a behavior fully corresponds with results published in the previous 
researches [8, 6, 9, 16]. Since the parameter /3 corresponds (in the thermodynamic model interpretation) 
to the inverse temperature v discussed in the section 2 we can now provide a deeper insight into changes 
of a vehicular microstructure. Contrary to the previous investigations, in this paper we present (for the 
first time, as we believe) a more-detailed study of time-clearance distribution that is performed not in 
density windows only (as for example in [16, 8]) but in flux-density windows. The presented analysis shows 
that the principal quantity determining a shape of headway distribution is the traffic density. Deviations 
caused by changes of a flux J (if density g is fixed) are almost imperceptible. Moreover, estimated values 
of a are low which confirms the results published formerly. If analyzing the minimum distance x(Q!, /3) we 
can see (figure 8) that the tested hypothesis is confirmed predominantly in the region of pure congested 
flows. A red triangle (boundary of a color area in figure 8) characterizing windows W{g,J) with the 















Inner structure of vehicular ensembles and random matrix theory 


11 


highest deviations can be explained as follows. The east and south sides of the triangle indicate regions 
where the number of data is not sufficient (and histograms are not therefore smooth enough), whereas 
the west side represents windows where congested traffic is mixed with free-flow traffic. 
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Figure 8. Color representation of a minimum statistical distance x{ct, fi) 
obtained by the MDE-method. We plot values x(“MDEiD ivide) in various 
flux-density windows W{i>,J). Squares symbolize locations of window-centroids 
{q + I^qI2, J + Aj/2). The plotted values are valid for a fast lane. 


For completeness, we note that the results pictured in figures 5, 6, 7, and 8 has been extracted from 
fast-lane data where car-car interactions are stronger (due to high speeds). Alternative results (valid for 
main-lane data) are visualized in figures 11, 12, and 13 in the Appendix 5.3. Comparing results obtained 
for fast and main lanes one can detect significant differences that are particularly noticeable in figures 6 
and 11. Here we detect different ranges of parameter /3 and also different maxima for flux and density. 
All these disparities are caused predominantly by lower speeds of cars moving in a main lane. 

Optimal parameters d, [3 have been estimated by the method MLE as well. However, their values are 
(similarly to the results described in the subsection 1.3) extremely close to dMDE,/3MDE and therefore we 
will not depict them in figures. 

3 . 4 . Damped matrices vs. freeway traffic - direct link 

Finally, we analyze a direct link between the random matrices investigated and traffic-flows monitored on 
expressways. For these intentions we minimize (analogously to the previous approaches) the statistical 
distance x calculated for a histogram of vehicular clearances and for a histogram of level spacings in 
DUE—ensembles. Thus, we search an optimal value of the damped parameter g for which that distance 
is minimum. Comparison of the above-referred histograms is illustratively visualized in the figure 5 (see 
histograms and crosses). Optimal values of g are plotted in figure 9. Associated minimum distance 
are plotted in figure 10. As understandable, the values visualized in figure 10 are higher that those 
enumerated as a distance between real-road distributions and the functional family (21). The main 
reason is that the former represents one-parameter fitting, whereas the latter represents an estimation 
that is two-parametric. 

4. Highlights, discussion and conclusion remarks 

This article has two scientific layers: mathematical (dealing with new types of random matrices) and 
physical (investigating properties of inner structure of vehicular streams). Firstly, this paper builds on 
the pioneering work [13] dealing with new classes of random matrices. It brings a substantially deeper 
view into a spectral structure of damped matrices than the article [13] itself. Specifically, using theoretical 
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Figure 9. Color representation of the estimated value g obtained by the MDE- 
method. We plot values x(^mdEi /^mde) in various flux-density windows W (p, J). 
Squares symbolize locations of window-centroids (^-l-A^/2, J-|-Aj/2). The plotted 
values are valid for a fast lane. 
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Figure 10. Color representation of a minimum statistical distance 
obtained by the MDE-method. We plot values x(“MDEi/^ mdb) in various 
flux-density windows W{(),J). Squares symbolize locations of window-centroids 
{q + ^qI2, J + Aj/2). The plotted values are valid for a fast lane. 
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approaches we predict an analytical formula for level density of DUE—matrices, describe particularities of 
the unfolding-procedure, and estimate specific values of parameters used in a predictive formula for the 
LS-distribution. Parameters in a two-parametric family (13) are estimated with help of the minimum 
distance estimator and the standard Kolmogorov estimator. 

Secondly, we show that the same mathematical formula may be used as a theoretically-substantiated 
prognosis for a time-clearance distribution describing statistics of time-gaps among consecutive vehicles 
in real-road traffic. Indeed, both the estimation-procedures used (minimum distance estimation and 
maximum likelihood estimation) prove that the function (21) represents a suitable candidate for a time- 
clearance distribution. In this work, for the first time, the time-clearance distribution is monitored in 
various flux-density sub-windows inside the entire flux-density plane. Thus, the output represents a 
map of estimated values a and P depending on traffic density and flux. Since the parameter P admits 
a thermodynamic interpretation (as is discussed in [7, 8]) the changes of the socio-physical coefficient 
P = P{g, J) reflect a level of mental pressure which the drivers are under during a given traffic constellation. 
As is evident in figure 6 for fast lanes and 11 for main lanes, a driver’s mental strain depends predominantly 
on actual density and only marginally on actual flow. Among other things, we can clearly detect here 
a process of traffic saturation accompanied by corresponding changes in a traffic-microstructure. More 
specifically, with increasing density the coefficient P is growing, which is reflected in more systematic 
configurations of vehicles (as a consequence of the fact that the variance of a clearance distribution 
drops). Another benefit for physics of traffic is a relatively small statistical distance (if compared with 
previously published analyzes) between real-road histograms and a theoretical prediction. It confirms a 
quality of the approach used. 

Finally, it has been demonstrated that the new type of random matrix ensembles ((/—parameterized 
class DUEg) may be used to describe statistical distributions of inter-vehicle gaps. To be specific, similarly 
to the Ref. [2, 3, 4] (where it is proven that time-gap distribution monitored for Cuernavaca buses is 
the same as the level spacing distribution of CUE- matrices) we now find that a similar correspondence 
exists between cars on expressways and the ensemble DUE. For a chosen traffic state (described by flux 
and density) there exists an optimal value of the parameter g which the LS-distribution corresponds for 
to time-clearance distribution measured on a expressway. Figure (9) shows that the force-parameter g 
is increasing with traffic density, which supports an intuitive notion on changes of interaction forces in 
thickening flows. 

To be sure that the above-mentioned similarities between traffic and Random Matrix Theory are 
more substantial (than presented in our article) one should show that the detected correspondence remains 
valid also if investigating second order characteristics (statistical rigidity or number variance). This will 
be covered in continuing researches. Preliminary numerical analysis, however, shows that a link between 
vehicular samples and eigenvalues of random matrices is definitely not accidental. 
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5. Appendix 

5.1. Convergence of the distribution (2) 

Let 10 be the space of all smooth functions with a bounded support (according to [20]) and 10' be the 
class of generalized functions (linear and continuous functionals over a space 10.) Then for an arbitrary 
ip{x) G 10 : 

r C(il) 

lim {q^)^e{x)]ip{x)) = lim / -©(e - |a;|) exp[—^- -]g}{x)dx = 

£-5-0+ e-s-0+ £ £^ — 

C(a9) lim / 0(1 - |j/|)exp[--- :^]ip{£y)dy = 

£^o+ l-y^ 

C{d)(f{0) [ exp[--^-^]dy = (p{0) = {S{x);ip{x)). (26) 

J-i 1 - r 
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Therefore lim£_).o+ q^^g{x) = 6{x). Analogously lim^^o+ for 


lira . — 

^^ 0 + \y'2TrK 


"^■^ipi^x) = lim 


^->■0+ Jr 


' 2 ^ ipix) dx = 


lim 

x^0+. 


1 


.-yV2 


ip{>cy)d.y = (^(0) = ((5(x); (/?(x)). 


/r 

For completeness, we add that ©(x) stands for the Heaviside step function. 


(27) 


5.2. Normalization and scaling in the distribution (13) 

^2 2 

— t.f.a normalization and scaling conditions for the density (13) 


Denoting J{x,>c) = /q°° e « e 
transform into the equation 


J{x, k) = — 


1 dJ 


(28) 


2m: Ok 

We add that the substitutions D = , (3 = x^/4 have been applied for convenience. Using the identity 

.j(x,M) = ^{^y+-ic Q-+i(x>c) one finds 

(2>c^ — 1 — a)lCa+i{xK) = —x>c/C(j_|_]^(x>c). (29) 

By means of the asymptotical expression 

2 z 


lCraiz)z^ (2m-2)!! ( 1 + 


2m — 1 


(2m > 1) (30) 


we obtain the scaling equation for small x : 

2x>f:^ + >c^(l + 2 q; — x^) — 2xM:{a + !) — (« + 1)(2q; + 1) = 0. (31) 

Near the infinity this equation can be simplified to the form 
4>c^ — 2xx — (2Qf + 3) = 0, 

which leads to a detection of the linear trend in the scaling constant 


D{P) ~ P + — + a (/3—>■ oo). 
On contrary, with help of the approximating expression 


ICmi.z) 


z—^oo / TT 


2 ® 


1 4m^ — 1 

—F= + 


(32) 

(33) 

(34) 


we obtain the scaling condition for large x : 

l6xM^+2M^{3+8a+4a‘^-4xy-XM{15+l6a+4:ay-i4:a^+12a‘^+lla+3) = 0.(35) 

Surprisingly, this formula leads (after a simplification) to the scaling formula (32) as well, and therefore 
to the same linear asymptotics (33). 

Above that, since the variant /3 = 0 in (13) represents the well-known gamma distribution it is easy to 
prove that lim,g_>.o+ D{p) = a + 1 and therefore the final scaling formula (phenomenologically corrected) 
reads (14). Finally, for fixed values a, /3, and D we have (15). 


5.3. Minimum distance estimations for main lanes 

In this part of Appendix we picture results of minimum distance estimations applied to traffic data 
extracted from main lanes of the Prague Ring. 
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Pmde 



Figure 11. Color representation of values (3 estimated by minimum distance 
estimation (23). We plot values /3 ]vide various flux-density windows W{q^J). 
Squares symbolize locations of window-centroids A^/2, J Aj/2). The 

rectangle at the bottom right schematically visualizes the window-size used. Gray 
dots correspond to locations {Qk'iJk) of individual data sub-samples before a 
separation procedure. Unstained areas indicate either a lack of data or belonging 
to a free-traffic regime. The plotted values are valid for a main lane. 
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Figure 12. Color representation of values a estimated by minimum distance 
estimation (23). We plot values log(d]viDE) in various flux-density windows 
W{q^ J). Squares symbolize locations of window-centroids (^ + A^/2, J -h Aj/2). 
The plotted values are valid for a main lane. 
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Figure 13. Color representation of a minimum statistical distance 
obtained by the MDE-method. We plot values x(diviDE?/^ mde) various 
flux-density windows W{q,J). Squares symbolize locations of window-centroids 
(q -h A^/2, J -h Aj/2). The plotted values are valid for a main lane. 
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